Orbital Ordering, New Phases, and Stripe Formation in Doped Layered Nickelates 
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Ground-state properties of layered nickelates are investigated based on the orbital-degenerate Hub- 
bard model coupled with lattice distortions, by using numerical techniques. The Neel state composed 
of spin S—l ions is confirmed in the undoped limit x=0. At £=1/2, novel antiferromagnetic states, 
called CE- and E-type phases, are found by increasing the Hund's coupling. (3a; 2 — r 2 /3j/ 2 — r 2 )-type 
orbital ordering is predicted to occur in a checkerboard-type charge-ordered state. At £=1/3, both 
Coulombic and phononic interactions are found to be important, since the former stabilizes the spin 
stripe, while the latter leads to the striped charge-order. 

PACS numbers: 75.30.Kz, 75.50.Ee, 71.10.Fd, 75.47.Lx 



The existence and origin of "striped" structures con- 
tinues attracting considerable attention in the research 
field of transition- metal oxides 0. In a system with 
dominant electron-electron repulsion, the Wigner-crystal 
state should be stabilized, but in real materials more 
complicated non-uniform charge structures have been 
found. In Nd-based lightly-doped cuprates, neutron scat- 
tering experiments revealed incommensurate spin struc- 
tures 0] where antiferromagnetic (AFM) spin stripes 
are periodically separated by domain walls of holes. In 
La2-rrSr a; Cu04, dynamical stripes are believed to exist 
along vertical or horizontal directions (Cu-0 bond direc- 
tion) |3(. In nickelates, the charge-ordered stripes are 
along the diagonal direction In manganites, evidence 
for striped charge-ordering also along the diagonal direc- 
tion has been reported in the AFM phase for x>l/2 
while short-range diagonal stripe correlations have been 
found in the ferromagnetic (FM) phase at x<l/2 6J. 

In general, stripes can be classified into metallic or 
insulating. In L^-zSr^CuOi, the dynamical stripes ex- 
hibit metallic properties, but they are easily pinned by 
lattice effects and impurities. In Lai.6- a ;Ndo.4Sr a; Gu04, 
stripes along the bond-direction are pinned by lattice 
distortions 1], but they are still metallic. Intuitively, 
vertical or horizontal stripes could be associated with 
the formation of "rivers of holes" , to prevent individual 
charges from fighting against the AFM background 0. 
Such stripes should be metallic, even if they are pinned, 
since they are induced by the optimization of hole motion 
between nearest-neighbor Cu-sites via oxygens. 

However, in the diagonal stripes observed in mangan- 
ites and nickelates, charges are basically localized, indi- 
cating that such insulating stripes are not determined 
just by the optimization of the hole motion. In the FM 
state of manganites, the hole movement is already opti- 
mal and, naively, charges should not form stripes. Ob- 
viously, an additional effective local potential must be 
acting to confine electrons into stripes. If such a po- 
tential originates in lattice distortions, it is expected to 
occur along the bond direction to avoid energy loss due to 



the conflict between neighboring lattice distortions shar- 
ing the same oxygens. Then, static stripes stabilized by 
lattice distortions tends to occur along the diagonal di- 
rection, as shown in the stripes of the FM-phase of man- 
ganites, stabilized by Jahn- Teller (JT) distortions ||. 

In simple terms, vertical or horizontal stripes in 
cuprates can be understood by the competition between 
Coulomb interaction and hole motion, while diagonal 
stripes are better explained as a consequence of a robust 
electron-lattice coupling. However, a difficulty has been 
found for theoretical studies of stripe formation in doped 
nickelates, since both Coulomb interaction and electron- 
lattice coupling appear to be important. Since the Ni 2+ 
ion has two electrons in the e g orbitals, on-site Coulomb 
interactions certainly play a crucial role to form spins 
5=1. When holes are doped, one electron is removed 
and another remains in the e g orbitals, indicating that 
the hole-doped site should become JT active. Then, in 
hole-doped nickelates both Coulombic and phononic in- 
teractions could be of relevance, a fact not considered in 
previous theoretical investigations. 

In this Letter, charge ordering in doped nickelates is 
investigated based on the orbital-degenerate Hubbard 
model coupled to lattice distortions, using numerical 
techniques. After confirming the Neel state composed 
of 5=1 spins at x=0, both cases x=l/2 and 1/3 will be 
analyzed. At x=l/2, novel AFM phases called CE- and 
E-type have been unveiled, which are consistent with ex- 
perimental results. Further including the JT-type coop- 
erative distortion, (3x 2 — r 2 /3y 2 ~ r 2 )-type orbital order- 
ing is predicted. For x—1/3, an incommensurate spin 
structure is induced by the Coulombic model, including 
the level splitting between e g orbitals, but the charge 
stripe does not appear. To reproduce simultaneously spin 
and charge stripes, it is important to include the strong 
coupling of e g electrons to lattice distortions, originating 
in the in-plane oxygen motions. 

The model for nickelates includes three important in- 
gredients: The kinetic motion of e g electrons, Coulomb 
interactions among e g electrons, and electron-lattice cou- 
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FIG. 1: (a) Spin correlation S(q) vs. A for x=0. (b) Two 
kinds of local e g -electron arrangements for x=0. (c) AFM 
spin pattern theoretically determined for A<3. 

plings between e g electrons and distortions of the NiOe 
octahedra. Note that the electron-lattice term is divided 
into couplings for the apical and in-plane oxygen motions. 
In layered nickelates, all NiOg octahedra are significantly 
elongated along the c-axis, splitting the e g orbitals. This 
splitting from apical oxygens should be included explic- 
itly from the start and, then, the in-plane motion should 
be studied. The Hamiltonian H is given by 

H = - * 77 '4 7CT d i+a 7 ^+ A H( nia ~ nib V 2 
ia77 'a i 

i,7 i^c' 

+ U> 5Z "ia^ib + J' 51 d i 7 T d i 7 l d i7'l rf i7'T: ( l ) 
i i,7^7' 

where di acr (diha) is the annihilation operator for an 
e g -electron with spin a in the d x 2_ y 2 (d 3z 2_ r 2) orbital 
at site i, n- na — d\ ia di~ ll7 , n^—^^n^u, a is the vec- 
tor connecting nearest-neighbor sites, and , is the 
nearest-neighbor hopping amplitude between 7- and 7'- 
orbitals along the a-direction, given by t* a =— V3t* h = 
— V3t* a =3t* h =3t / 4 for the x-direction and t^ a —V3t^ h = 
\/3iba— 3^=3^/4 for the y-direction 0. Hereafter, t 
is the energy unit. In the second term, A(>0) is the 
level splitting between a- and b-orbitals. In the Coulomb 
interaction terms, U (U 1 ) is the intra-orbital (inter- 
orbital) Coulomb interaction, J is the inter-orbital ex- 
change interaction, and J' is the pair-hopping amplitude 
between different orbitals. Due to the relations J=J' and 
U=U' '+J+J' the independent parameters are U' and J, 
with U'>J |9(. The calculations are all carried out using 
standard exact-diagonalization techniques. It is crucial 
to use this kind of unbiased methods for this first study 
that includes both Coulomb and lattice effects, even if 
the technique restricts us to small TV-site clusters. 

First, consider the undoped case. The calculation is 
done for an 8-site tilted cluster, equivalent in complexity 
to a 16-site lattice for the single-band Hubbard model. 
Since at all sites the two orbitals are occupied due to 



the Hund's rule coupling, the JT distortions are not 
active and it is possible to grasp the essential ground- 
state properties using H. In Fig. 1(a), the Fourier 
transform of spin correlations is shown vs. A, where 

with s?= E 7 (<t^t 

— dj 7 |C?i 7 jJ/2. As expected, a robust (7r,7r) peak can 
be observed for A<3, suggesting that the AFM phase 
is stabilized by super-exchange interactions. The rapid 
decrease of >S(7r, tt) for A>3 is understood by compar- 
ing the energies for local triplet and singlet states, as 
shown in Fig. 1(b). The ground-state properties change 
at U'—J=U—A, leading to A=3 J for the transition. The 
spin structure at a;=0 is schematically shown in Fig. 1(c). 

Let us turn our attention to the case x=l/2. The 8- 
site tilted lattice is again used for the analysis, and the 
phase diagram Fig. 2(a) is obtained for A=0.5 10]. In- 
creasing J, an interesting transformation from AFM to 
FM phases is found. This is natural, since at large J 
the system has a formal similarity with manganite mod- 
els, where kinetic-energy gains lead to ferromagnetism, 
while at small J the magnetic energy dominates. How- 
ever, between the G-type AFM for J«0 and FM phase 
for Jw[/', unexpected states appear which are mixtures 
of FM and AFM phases, due to the competition be- 
tween kinetic and magnetic energies. Typical spin corre- 
lations S(q) are shown in Fig. 2(b). Note that peaks at 
q=(n, 0) and (7r/2,7r/2) indicate "C" and "E" type spin- 
structures, respectively (the notation is borrowed from 
the Mn-oxide context |9|). Double peaks at q=(7r,0) and 
(7r/2,7r/2) denote the CE-type structure, frequently ob- 
served in half-doped manganites 01 ■ ^ n half-doped nick- 
elates, the CE-phase is expressed as a mixture of type (I) 
and (II) in Fig. 2(c), depending on the positions of the 
S=l and S'=l/2 sites, although the "zigzag" FM chain 
structure is common for both types. The E-type phase 
is also depicted in Fig. 2(c). Note that the charge cor- 
relation always exhibits a peak at q=(7r,7r) (not shown 
here), indicating the checkerboard-type charge ordering. 

In experimental results, a peak at (tt/2,tt/2) in S(q) 
has been reported 0, suggesting an AFM pair of 5=1 
spins across the singly-occupied sites with holes. More- 
over, the checkerboard-type charge ordering has been ex- 
perimentally observed Thus, the spin-charge patterns 
of CE(II)- and E-type are consistent with the experimen- 
tal results. Our phase diagram has a robust region with 
a peak at (tt/2,tt/2), both for CE- and E-type phases, 
although the CE-phase exhibits an extra peak at (ir, 0). 
Whether the E- or CE-phases are present in nickelates 
can be studied experimentally in the future by searching 
for this (71", 0) peak. Note that if diffuse scattering exper- 
iments detect the AF correlation along the hole stripe, 
as has been found at x=l/3 12], the CE(II)-type may be 
the only possibility. Summarizing, the spin-charge struc- 
ture in 2;= 1/2 experiments can be understood within the 
Hamiltonian H by assuming a relatively large J. 
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Consider now the effect of in-plane oxygen motion (api- 
cal oxygen motions have already been included as an e g - 
level splitting). Assuming that oxygens move along the 
Ni-0 bond direction, the extra electron-phonon coupling 
term is written as 

Heph = g y^[-Qli(Tlia + mb) + Q-ZiTxi + QsiT^i] 
i 

+ (k/2)J2(PQli + Ql + Qti), (2) 

i 

where g is the electron-lattice coupling constant, Qu is 
the breathing-mode distortion, Q 2 i and Q^i are, respec- 
tively, the (x 2 — y 2 )- and (3z 2 — r 2 )-type JT distortions, 

T *i = T,a( d L* d ibv + 4b* d i^)> and T zi = T,a( d La d i^ ~ 

d iha diba)- The second term is the quadratic potential 
for adiabatic distortions, where k is the spring constant 
for the JT-mode and (3 is the spring-constant ratio for 
breathing- and JT- modes. From our experience in man- 
ganites, this ratio is here fixed to j3=2 |13|. 

Since all oxygens are shared by adjacent NiC"6 octa- 
hedra, the distortions are not independent. To consider 
such cooperative effect, in principle, the O-ion displace- 
ments should be optimized. However, in practice it is 
not feasible to perform both the Lanczos diagonalization 
and the optimization of all oxygen positions for 6- and 
8-site clusters. In the actual calculations, Qii, Q?n and 
Qsi are expressed by a single parameter d, for the shift 
of the O-ion coordinate. Note that the unit of d is g/k, 
typically 0.1~0.3A. Then, the total energy is evaluated 
as a function of d to find the minimum energy state. Re- 
peating these calculations for several distortion patterns, 
it is possible to deduce the optimal state. 

After several trials, the optimal distortion at x—1/2 is 
shown in Fig. 2(d). The diagonalization has been per- 
formed at several values of d on the 8-site distorted lat- 
tice and the minimum in the total energy is found at 
d=0.3 (Fig. 2(e)). Here, the dimensionless coupling con- 
stant A is defined as X—g/Vkt. As mentioned above, 
even without -ff cp h, the checkerboard- type charge order- 
ing has been obtained, but the peak at g=(7r,7r) signifi- 
cantly grows due to the effect of lattice distortions. Note 
that the distortion pattern in Fig. 2(d) is essentially the 
same as that for half-doped manganites. This is quite 
natural, since JT active and inactive ions exist bipartitely 
also for half-doped nickelates. Then, due to this JT-type 
distortion orbital ordering for half-doped nickelates is pre- 
dicted, as schematically shown in Fig. 2(d). The shapes 
of orbitals are determined from the orbital densities, (r z i) 
and (Txi) (Fig. 2(f)). The well-known alternate pattern 
of 3x 2 —r 2 and 3y 2 —r 2 orbitals in half-doped manganites 
is denoted by dashed lines. Increasing d, the shape of 
orbitals deviates from 3x 2 —r 2 and 3y 2 —r 2 , but it is still 
characterized by the orbitals elongating along the x- and 
y-directions (see insets of Fig. 2(f)). It would be very 
interesting to search for orbital ordering in half-doped 
nickelates, using the resonant X-ray scattering technique. 
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FIG. 2: (a) Ground-state phase diagram at x—1/2. (b) S(q) 
for the CE- and E-type phases, at the couplings indicated, 
(c) Spin and charge patterns for the CE- and E-type phases. 
These are schematic views, since local charge-densities in 
practice are not exactly 1 and 2. (d) Numerically obtained 
cooperative distortion pattern for an 8-site lattice at a;=l/2. 
Black and open circles indicate Ni and O ions, respectively. 
Open symbols indicate e g orbitals in the optimized state, (e) 
Total ground-state energy vs. d for a;=f/2. (f) Orbital den- 
sities (r z i) and (r x i) for sites 1-4. See (d) for the site labels. 
Optimized orbitals at d=0.3 for sites 1 and 3 are also shown. 



Now let us move to the case cc=l/3. If the actual 
expected stripe structure at x=l/3 is faithfully consid- 
ered 0, it is necessary to analyze, at least, a 6x6 cluster. 
However, such a large-size cluster with orbital degeneracy 
cannot be treated exactly due to the exponential growth 
of the Hilbert space with cluster size. Then, a covering 
of the two-dimensional (2D) lattice using zigzag 6-sites 
clusters is considered (Fig. 3(a)) by assuming a periodic 
structure along the diagonal direction. The phase dia- 
gram obtained by analyzing the zigzag 6-site cluster for 
H is in Fig. 3(b). Typical spin and charge correlations 
are in Figs. 3(c) and (d), where C(q)=(l/N)J2 iti e i? ' (i ^ 
((m - (n))-(nj - (n))), with n ; =^ 7 n i7- 

Since the momentum q is defined along the zigzag di- 
rection, the phase labelled by q=2n/3 in Fig. 3(b) denotes 
an incommensurate AFM phase with the proper spin 
stripe structure. The phase labelled by q—n/3 indicates 
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FIG. 3: (a) Zigzag 6-sites cluster covering the 2D lattice. 
Black circles denote Ni ions, and dashed lines indicate hole 
positions, (b) Phase diagram at x=l/3. Each phase is char- 
acterized by the momentum that shows a peak in S(q). (c) 
S(q) and (d) C(q) vs. J for U'=6 and A=0.5. (e) Coopera- 
tive distortion pattern for the zigzag 6-sites cluster at 1=1/3. 
(f) Total ground-state energy and (g) C(q) vs. d for x=l/3. 



a spin spiral state, which will eventually turn to the FM 
phase in the thermodynamic limit. Thus, the spin stripe 
phase appears between the commensurate AFM and FM- 
like phases, similar to the case of £=1/2. However, as 
seen in Fig. 3(d), C(q) in the spin stripe phase does not 
show the striped charge structure (q=2n/3). Rather, bi- 
partite charge ordering characterized by a peak at q=ir 
still remains. Namely, the Hamiltonian H can explain 
the spin stripe, but does not reproduce the striped charge 
ordering at £=1/3, indicating the importance of iZ C ph- 

Consider now the effect of i? e ph for £=1/3. After eval- 
uating total ground-state energies for several kinds of 
distortions, the pattern in Fig. 3(e) has been found to 
provide the optimal state at £=1/3. This type of distor- 
tion induces a spatial modulation of the level splitting as 
— <5i/2=r5 2 =f5 3 =— 64/2=65=6^ Q, where Si is the level 
splitting caused by the in-plane oxygen motions, and the 
site numbers are in Fig. 3(e). The minimum energy is 
found at d=0.1 (Fig. 3(f)). The modulation of level split- 
ting stabilizes the striped charge ordering characterized 
by a q=2n/3 peak in C(q) (Fig. 3(g)). 

Note that (3x 2 — r 2 /3y 2 — r 2 )-type orbital ordering does 
not occur in Fig. 3(e). Phenomenologically, such or- 



bital ordering tends to appear in a hole pair separated 
by one site, the unit of the "bi-stripe" of manganites 
.j5j. However, such a bi-stripe-type ordering contradicts 
the £=1/3 striped charge-ordering, and the bi-stripe-type 
solution was found to be unstable in these calculations. 
One may consider other distortion patterns which satisfy 
both (3x 2 — r 2 /3y 2 — r 2 )-type orbital and striped charge- 
ordering, but in such distortions no energy minimum was 
obtained for d>0. After several trials, Fig. 3(e) has pro- 
vided the most optimal state. 

Summarizing, possible spin, charge, and orbital struc- 
tures of layered nickelates have been discussed based 
on the e g -orbital degenerate Hubbard model coupled 
with lattice distortions. To understand the nickc- 
late stripes, both Hund's rule interaction and electron- 
lattice coupling appear essentially important. At £=1/2, 
(3a; 2 — r 2 /3y 2 — r 2 )-type orbital ordering similar to that 
in half-doped manganites is predicted. Even FM phases 
could be stabilized by chemically altering the carrier's 
bandwidth. For £=1/3, a spatial modulation in level 
splitting plays an important role for stripe formation. 
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